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INTRODUCTION 


The theory of LSC is а special edition of the theories of adsorp- 
tion from multicomponent-liquid mixtures on solid surfaces, when опе 
ot the components of the mixture (solute) appears іп a very small 


concentration, This adsorption process is governed by several physical 
factors, the most important of them being: 

1, Тре interactions between solvent and solute molecules in the 
surface and the ^ulk phase; 

2, The differences in the surface areas occupied Бу different 
molecules; 

3. The effecte of multilayer adsorption; 

& The energetic heterogeneity of solid surfaces, 
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There is no generai theoretical solution for this problem, Howates, 
several attempts have been made: most of them based on lattice theories 
of solutions. Surprisingly, their application to the specific case of LSC 
was very small L1-3Z. There may be several reasons for that: 

Firstly, the simple theories of LSC C4-18] were good enough 
to describe the retention data at the common degree of the experimen- 
tal accuracy. Secondly, the more advanced theories of LSC must still 
be tractable for the majority of the scientists working in LSC, 

This publication presents the theory of LSC, which takes into 
account all the basic physical factors enumerated above, This theory, 
is applied successfully to describe quantitatively the retention data 
т the typical LSC systems, investigated by Siaats et al. C131. 


THEORY 


Let us consider firstly the competitive adsorption of the molecules 
A and B in the i-th lattice plane, Assuming that the molecules occupy 
the same surface (lattice) area, this competitive adsorption can be 


expressed as the following quasi-chemica! reaction, 


gj 3 (1) , 
(в) (6) 1 (д) (9) ша (в) (9). a) (2) 
The related condition for the thermodynamic equilibrium reads, 
(ь). ит Wy r6) ; 
^g "id ЗА (2) 


where [A (1) апа © ) аге appropriate chemical potentials in the 
bulk phase, and in the ith lattice Plane, Denoting by the symbols x 
and (Y appropriate mole fractions: and activity coefficients and using 


the relation 


(== Эл + кт x (3) 


one can rewrite equ, (2) to the following equivalent form, 


а) роб) „ү ©) 


— Bi = K (4) 
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where 


рые о AO oup, 
PNE js kT кт - У 


Equation (3) can be rewritten to the foilowing Langmuir-like form 


C. L L wt ELENA L с 
4 (Ya ДУХ * Kear "B Qs IN вд 
‘Now, let us consider the case, when the molecules A and B 


occupy different surface areas c, and e pr Then, eqs. (1-5) take 
the following form, ` 


(6) 


TP 28 фа w, et f? n (7) 

D а, е0), SA pg lO с 
(1) (1) (e) (5) Y мА 

"a E ї a e Коко (9) 
x хв $ в 
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Let us consider now the simultaneous competitive adsorption of 
а third component C, which we shall assume to be the solute appear- 
ing in very small concentrations, Rs competitive adsorption in. the i-th 
lattice plane will be described by the following set of the quasi-chemi- 
cal reactions: 


(7 , ас) S (ay T Ace » (a1) 


^2 


if, ze (cf) S (wl, Ee . tiz) 
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When “д = c в = с: the statistical thermodynamics yields the 


following equation for the adsorption isotherm of C in the presence of 


А and в 

1 aem o» 
хер 0D. РУ (yy а) 7} GY €)' 
En Ac. (fc {сл E КАСА, (| А. ТУНХ Kuci*p. wd Qo) 


Moreover, since х0) is very small, the above equation reduces to 


the simpler form, 


«Du p) 
L 


I Ы (14) 
х l ^ 
BCJ"B' - B ви 


АСИХА x і AJ UM 
Till now, there is no general solution for the case when A at 
е в Я о с’ bes for adsorption isotherm from multicomponent (ternary) 
liquid mixtures of molecules occupying different surface (lattice) areas. 
Below, we are going to proposa some approximate but general solution 
© of this problem, Let us remark to this purpose that the denominator 
in equ, (14) has the following meaning, 


E PEE sc ENT EE 


Ге {ДУ ау (1) х (y OO 
| А А хв в А dip 
ie ыргуу, a ter e e Eic tay 28) 
Ai Bi ХА Sai 


Because of its small prr уте Qr) 29) of solvent C will 

; L)y (L €) у (€) 
not affect much the ratio ( XA ФА Ix si Ad 
assume that the deviations of ( aoc) or (ese) from unity 


could be taken into account in a form of a correction in equ, (14), 


). We may, therefore, 


which will remain basically unchanged, Writing the equations for the 
thermodynamic equilibria related to the quasi-chemical reactions (11-12) 


seems to suggest, how to take these corrections into account: 


(1), *c (9) (c) “с. (x) 
Me а МОВА. ke и? , (16) 
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We can see that the deviations СИР ) ana (< Че) 
from unity cause the solvents A and В to exhibit some apparent 
effective value of their chemical potentials in the quasi-chemical 
reactions, in Бом bulk and the surface phase, These apparent values 
of ex and (^^ (9) аге equal to their true values, multiplied by 

= 
the factors( af A A) and (Ilp) 

Now, we are coming to the essential point of our approximation, 

Namely, we shall keep the basic form of ee Sm ой the terms in 
(u сууз 
the denominator wil be г ye by к $ A I E v A and 


ДК 


Doing so, we arrivo E У followings form of the ‘distribution coefficient 


Kee . of solute, 


> according to whet we have” аа above, 


Сы | | 
E | (өрү) Ga 
E өү 55 ЕВ 
ка >й? Oy ©) А x1 X Gu EDI 
The experimentally measured distribution coefficient Mes represents 


some average over ай the adsorbed layers (lattice planes), where 
the concentration of С is different from xc" „ However, because of 
the smali concentration of C, and the usual condition that Won » 
pun Kc the solute C will practically be adsorbed only 
in the first layer. At the same time, however, the Solvents may well 
be adsorbed in a multilayer fashion, Such a model of surface phase 
will be accepted in our further consideration, а 

Let us assume. юг the moment that the distribution of the solute 
C between the surface, and the bulk phase is caused only by the 
interactions between the solute and the solvents molecules, In other 
words, we are staying purely on the ground of the “interaction” model 
of HPLC, The difference between the concentrations of the solute in 
the surface, and the bulk phase is caused by the different concentra- 
tions of the {A+B) solutions in the surface and tne bulk phase, in 
which the solute is "soluted", In such a refined form, the "interaction" 
model] has been formulated firstly by OSci k L41 The distribution 


coefficient ae ee for this model is defined as follows 
* 
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(6) (1) 
x яс c 2 (19) 
cC C MEN y (<) 
^g C,i 
Thus, in our general treatment, the distribution coefficient TK 1 
2 
can be written in the following form 
x x се. $ > 
Koa “4 Keun Ke com * ; (20) 
where ХК. is the "competitive" distribution coefficient, 
C,com | 
SG “с 
— — m -1 
СА po “в 
"к - Ke x (1) S * ka MI B . {21) 
C;com |^АСл: ГА { ©) вс,1 ^B I (©). 
4,1 B,1 
In the absence of intermolecular interactions, Бе. when п, 1, 1. 


х the total distribution coefficient ca reduces to the "сотре С" 
factor MG om Beginning from now, we shel drop the subscript 
"1, іп ЖОЛА since we assume {пе monolayer adsorption of solute С, 

et us consider now another important physical factor which 
affects the distribution of solute between the surface, and the buik 
phase. This is the energetic heterogeneity of real solid surfaces, 
which causes the adsorption sites to be distributed among various 
adsorption energies. In our notation, it means some dispersion of 
the values of the parameters Kacy and Каса” ^ 

We shall assume further a random model of surface topography 
ie. the lack of any spatial correlationa between adsorption sites 
Ac, а Kic r Let X А(КАС, 1) 
represent the differential distributions of adsorption 


exhibiting the same value of К * 
and x gU. e 

Я SME K x i 
sites among various values of K AC. end кве, 1 Assuming that 


these functions are normalized to unity, we obtain 


^c 
x а a t bj *2A 
C,com dK c ] нс: баса КА VO * 
x ed 


j - 
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“с 
sf в"? 

a (agn 15) А) 
БеЗ Mal Kac) Лв(Квс,1 (22) 
oe Qa and B denote the range of variations of АВЕ and 
S 86,1 


The hitherto investigation of the competitive adsorption from 
binary liquid mixtures on solid surfaces, have shown, that tne 
distribution function X is in general а gaussian-like function, We 
shall approximate it further by the following analytical function, 


o 
1 а ied s ) 
Joe m 


> h wd (asa (23) 


where the parameter r describes the width of the distribution X (y) 
function, whereas y? is tne most probable value of the variable y. 
In the limit г — 0, AG) degenerates into Dirac delta distribution 

J (y - y?) , related to an ideally homogeneous surface, characterize 
| ed by: у = y? . The analytical form of the distribution function (23) 
makes the following expansion possible, for a function Е(у), 
continuous with all its derivatives, 


409 * Ч 

eae 2n . 
| Р(У(У) dv» к(у°)+г) BS о-в, [oe (24) 
а nel ayen y? 


where B, is Bernoulli's number, The investigation of adsorption from 
binaries on solid surfaces have shown, that the distribution function 
is, in general, a narrow one, That means, that r can, in general, be 
treated as a small parameter, and expansion (24) can be cut after 


the second quadratic term, Doing so, we obtain, 


+ ° 
22 2 
| PO) Доду em +7 i E (25) 
y o 


et Уу 
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Let us remark yet, that when K (y) is narrow, one can safely 
replace the integration interval (- Ф, + 29) by a finite domain $2 
such that y? € Gà. Further, according to the random model of surface 
topography, eccepted here by us, we can apply the expansion (24) 
when evaluating the multiple integral (22). Then, neglecting in this 
expansion the terms higher than quadratic with respect to the para- 


meters r's, we arrive at the following equation for хк 


C,com 
2 2 
x X30 Л (х o 
Кс, сот C,com + 3 KC com : 
к (26) 
ыле E 
үс) ‚^ (Od “в 
2 (x ta) ‚а (0) в ) ; 
АЕ А © Bc\*B [3 
А,1 в,1 
In equ, (26), ie com Means the competitive distribution coefficient 
LI 
for an ideally homogeneous surface, 
A a 
y (1) x) т 
хо Фо (1) "A Ф A 4:0. СЫ) Б В |-1 (27) 
C,com MC, д Г (6) BC "B у $ 
A, B,1 
whereas r and r are the parameters, describing the dispersion 
AC ed Wo Kò 
Li s 


A a 
of the coefficients КАс,1 $ Квс,1 ‚ апа К ACA 5 Кас 1 are their 
most probable values, 

In a similar way, one can take into account the effects of surface 


heterogeneity т the competitive adsorption of the solvents A and B: 


22 2. «S 
EE e AR DO 
е t TG | L (28) 
| 90g) NC ac e 


BA,1" “BA,1 


В га luated f (9), taking к“ KS d 
where хдт is evaluated from еди, ^ ng Кд = "TER 


ТАЗ is appropriate heterogeneity parameter, 

Now, let-us consider at last the problem of activity coefficients 
in both the adsorbed and the bulk phase. т the case of bulk activity 
coefficients the problem is relatively simple: Here, the power expansions 


for the logaritim of the activity coefficients are applied 
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n Hn. antag)? + v (xfi?) с Au, (29) 


in yi. а NOU * bg i x)? А с(ж00))* Use ’ (30) 


The coefficients а, b, c, are related to the molar excess free enthalpy 
of mixing aem) » by various theories of solutions, Let us consider, 
for instance, the very popular Redlich-Kister theory, which gives the 


following expression for Ban b " 


; . 2 
К owe x 0, () А + в(х С). x0). с(х (2). x) 
(31) 
3 
* p( (9?) - х) с ae -] . 
Application of the Rowlinson's relationship to equ, (31), yields tne 


following relations between the coefficients a, b, c and A, B, C: 


A + 3B + 5C + 7D a, = А = 3B + 5C - 7D 


у Чура в 
ba = -4(B + 4C + 9D) bg = 4(B = 4C + 9D) 
(32) 
сд = 12(C + 5D) св = 12(с - 5D) 
da = «32D da = 32D , 


Thus, knowing the activity coefficient of one component of the 
mixture (A+B), makes it possible to calculate easily the other one, 
by solving the linear system of equations (32). 


In the case of the solute C, we approximated its bulk activity 
coefficient 1n y by the following equation, 


m Fore Ac x + Agar + Beal 000)? ees x O2) 


вах; (№) 
ЗАВ. (33) 


in which the quadratic terms have to account for the deviations from 


a regular solution behaviour, as described by Scatchard-Hildebrand 


theory. However, for aeg) we used the quite general expression, - 
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see) бү, 0 а р). (34) 


We shail show it in ‘the next section, that eqs, (33-34) offer а very 
' effective approximation for In $ e А | 

Now, we are coming to the very essential] point of. our conside- 
ration, concerning the activity coefficients of the solvents A and B, 
and the solute C in the adsorbed phase, Everett suggests [19], 
that tne deviations from ideality in the adsorbed phase should. be 
fairly comparable with those in the equilibrium bulk phase, Following 
this suggestion, we approximated in eU . т. $ i and In is 
by the same equations (29), 30), (33) and (34), except that bulk 
mole fractions x » were replaced there by some locally averaged 


surface concentrations E: 


e » рхү°) кп ЕЗ * ч | Р ab 


‘ 


: In equ, (35) p is the fraction of the nearest neighbours lattice 
(adsorption) sites in the same lattice plane, whereas m is that 
fraction in the lattice planes below or above the considered laitice 
plane, In the case of the first lattice plane, i€. when i=1, it is 
assumed that the lattice plane lying below is empty. 


NUMERICAL RESULTS AND DISCUSSION 


For the purpose of iliustration we will consider the experimental 
data by Slaats et al, [13] we analyzed these data numerically 
,neglecting the effects of the non-ideality of the adsorbed phase, as 
well as the effects of the surface heterogeneity. With these simplifi- 
cetions we were not able to describe well the capacity coefficients 
of the solutes: Nitrobenzene and benzyl acetate, chromatocraphed in 
the two mixed mobile phases: (n-heplane 4 2-propanol) and (n-heptane 
+ ethyl acetate), on Partisil 10, Below, we are going to prove, what 
kind of agreement between theory and experiment can be obtained on 
the ground of our present theory, taking the effects of the nor-ideality 
in the adsorbed phase, and the effects of surface hetereogeneity into 
account, 
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; we bna with the activity coefficients of the solvents and . 
м in the bulk phase, which were megsured independently by 
Sleats et al [1з]. Our numerical investigation has shown, that in 
the case of the mixture (n-heptane + 2-propanol), the activity 
coefficient of 2-propanol| may well be represented by the polynomial 
‘of the seventh order, 


in Yi? - x22 (4 (17 13e (x 09)" А 2,92 (x 023) + 
5 6 ? (36) 
A 0.196 (x (2) + 0,628 (x +) + 8.0(х (9) 


where xí! ) is used to denote the mole fraction of n-he plane, in the 


case of die solvent у (n-heptane + ethyl acetate) a satisfactory © 


approximation for 1л т эк, can be obtained Бу the following fifth-order 


: 3 4 5 
in ROREM -3.4 (x (99) «4 A7(x 7) 2 (37) 


Figures. 1 and 2 show graphically ine effectiveness of the approxima- 
tions (365) and (37). In the next Figures 3 and 4, the effectiveness 

ot the approximation - is demonstrated, which can be obtained on the 

ground of оог eqs, (33), (34), for the Бы activity coefficients of the 
solutes, Table 1 reports the related coetticients Aca Ас Beas 
Bog found by fitting best eqs. (33-34) to te experimentel data by 
Siaats et aL [13]. 


Table 1 
Solute f ACA Асв Bea Bes 
(n-heptane + 2-propanol) 
Nitrobenzene 1.05 1,54 2.54 0.358 
Benzyl acetate 2.81 2,14 -0.349 0,325 
{n-heptane + ethyl acetate ) 
Nitre benzene 0816 0,932 -2,39 1.5 


Benzyi acetate -1.51 1.22 0,635 0,706 
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Now, let us consider at last the essential problem of the 
effectiveness of the approximation for the experimental capacity 
- coefficient kc» which can be obtained on the ground of our theory 
developed in the previous section, This experimental quantity has 
the following relation to the distribution coefficient "Kos 


vis) x 
кс = үт) Koa (38) 
where vis) and vim) are the volumes of the surface, and the mobile 
phase respectively. Р 

The theory developed Бу us introduces a number of parameters 
which, theoretically, could be measured in an independent way, 
especially designed experiment, For instance, all the parameters 
K 5A, 1* РАВ • (о „/ a) describing the competitive adsorption of 
solvents can be found by appropriate numerical analysis of the 
excess adsorption isotherm of solvents, measured by one of Те. 
Static methods [20], However, even in such a complete chromato- 
graphic experiment, as that carried out by Slaats et al, C13]tnis 
independent measurement of the excess isotherm of sotvents was 
not performed, j 

The next necessary step in our theoretical treatment is to 
decide which parameters play the most crucial role, and which ones 
are of a smaller importance. This has the purpose to decrease the 
number of ihe parameters, to be found by fitting best our equations 
to experimental data, 

Our extensive numerical investigation has brought us to tne 
conclusion that the parameters describing the competition between 
the solute and the solvents on the surface, are the most essential 
ones, И should, however, be noted, that no essential improvement 
is achieved when.one distinguishes between Fac and Tac * compared 
to Situation when they are assumed to be equal, As to the competitive 
,adsorption of solvents, only the parameter Kean affects strongly 
the agreement between the theory and experiment, The best-fit calcu- 
lation has appeared to be insensitive to the parameters РАВ апа 
( & leg). Of course, also the interaction parameters may оіау 
a more эг jess important role, but we assume that we know them 
from an independent bulk experiment, 

Looking tor & possibility to decrease the number of the best 
“АИ parameters, and to ‘the results of our numerical investigation, 
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we propose to use (at least for the. chromatographic Systems under 
investigation), the following simplified version of our theory, obtained 


with the assumption that “А = Ap” Ao =r -» г, and 


FAC BC 


rap = 0. Then, the basic equation (20) of our theory can be written 


in the following form 


2 AS 
z - 
" (d?) кдд AO MES TEM (1) a) 
cT Р (сб) Po КВА = Сад yo 2i 
Ca А,1 (38) 
x ся 
4) 
MC жа SUM 
B Хх &) 
B, 
where 
а (т) 
Pc F^ hey o Xn (39) 
yi) 
г : (40) 
iia 
| Ас. | dg 
Кл = (5400 [4 a EAM 
BA Toy 


{ Y A | 
A T = + Kp i(*5 y) 2 ] р (41) 
А1 | | 7 


In this мау, we have arrived at а simplified version of our 
theory which includes now the following five parameters to be fitted 
best numerically: PC, m, Kpa’ ( dcla) and г’. 

: Now let us explain the strategy of our best-fit calculation: 
Namely, for a given pair of solvents (A+B), we fitted simultaneously 
equ. (39) to the experimental data for both Nitrobenzene and Benzyl 


acetate, assuming that the parameter K must be the same, 


BA,1 
whereas the other parameters may be different, The results of this 
best-fit calculation are shown in Table 2, à 

Further, since we assumed that the competitive adsorption of 


solvents may have a multilayer character, we had to evaluate the 
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„concentration profile of the Solvents near the surface. To do it, we 
assumed that the surface forces of the adsorbent do not extend 
over the first adsorbed layer, ie, that for the second and higher 
layers K 3A = 1. Having the mole fraction of the solvents in an i 
adsorbed layer, say i, evaluated, we couid find the sae a 
excess rites of the more active ѕбіуепі B, 

A = xf) - ! (42) 
The total adsorption excess пе) which could be measured by using 
static methods, is just the sum of the excesses пр “= over all 
the adsorbed layers divided by the number of layers, The evaluation 
of the concentration profile had to be done in every step of our best- 
-fit calculation, Figures 5 and 6 show the concentration profiles 

(adsorption excesses) corresponding to the best-fit parameters 
collected in Table 2, | 


Table 2 
Solute P | т k^ d [o r^ 
i à Ў с BA,’ °С 
(n-heptane + 2-propanol 
Nitrobenzene 1.97 0.416 300.0 0.536 0,649 
Benzyl acetate ey E 0,416 300,0 0,402 . 0.248 
(n-heptane + éthyl acetate) 
Nitrobenzene _ 0.926 0.500 . 317.0 0.236 0,106 


Benzyi acetate 0,422 0,500 317.0 0,418 0.229 


` 


It can be concluded from these figures that the competitive adsorp- 
tion of the solvent mixture (n-heptane + 2-propanol) has a strong 
multilayer ‘character, whereas in the case of the mixture (n-heptane + 
+ ethyl acetate) a double-layer model should be quite satisfactory, 

High values of the parameter m mean that the vertical interactions in 
the adsorbed phase are predominant, At least in the first adsorbed 
layer which: makes the largest contribution to the total surface excess. 


Two things may be responsible for that: Firstly, ihe adsarmtion centers 
x К " 
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on Acrosil surface may lie at distances, at which two solvent molecules 
adsorbed on two neighbouring centers do not exhibit more some stronger 
interactions. Secondly, there may exist some orientation effects induced 
by the surface, which make the vertical interaction preferential, 

The values of C4 X4 ), which deviate much from unity suggest 
strong orientation effects of the solute molecules in the adsorbed phase. 

Figures 7, B, 9 show the agreement between experimental capacity 
coefficients measured by Slaats et aL [13], and our, theoretical ones, 
évaluated with the parameters presented in Table 2, The agreement 
seems to be quite satisfactory, ў А 

Figure 10 shows separately the contribution ke to the capacity 
coefficient ko’ due to the effects of surface heterogeneity 


L) к 


d ( ^ л БР" 
kc =: CM p Ка. 
EC e | ' 
, (43) 
2M a Ss 
{ (L) e { (1) a 
TCO Беде). 
А,1 B,1 


for the case of Nitrobenzene, Similar pictures are obtained for other 
systems, It сап be concluded from these figures, that the effects of 
Surface heterogeneity are especially important at very small concen- 


trations of the more active solvent 


CONCLUSIONS 


In this paper we have presented the theory of retention in LSC, 
based on a lattice model of the surface phase, In our theoretical treat 
ment, the expression for the capacity coefficient is written as a product 
of two terms: one of them corresponding to what was called the 
"interaction" theory in LSC, and the other one corresponding to the 
"competition" theory in LSC, Ң should, however, be noted that the 
"competitive" factor is also influenced by the interaction effects т 


‘both the surface ага һе bulk phase, On the other hand, the 
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"interection" term depends, (through Surface activity coefficients), 
on the mechanism of the competitive adsorption on the surface, 

Our numerical investigation has shown, that a satisfactory 
‘description of retention data in LSC requires some few basic charac- 
teristics of a LSC system to be taken into account: 

1. The interactions in the surface and the bulk phase, 

2. The differences in the surface areas occupied by solvent 
and solute molecules; 

3. The multilayer adsorption of solvents; 

4. The effects of the energetic heterogeneity of the solid 
(support) surfaces, 

In some cases, however, several simplifying assumptions can 

be accepted, They depend on the features of the particular LSC 
system under consideration, 
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STRESZCZENIE 


Myprowadzono réwnanie opisujące zależność współczynnika po- 
dziatu w adsorpcyjnej chromatografii cieczowej od składu binarnej fazy 
ruchomej, Uwzględniono nieidealność fazy ruchomej, różnice w rozmiarach 
cząsteczek substancji chromatografowanej i rozpuszczalników, energetycz- 
па niejednorodność powierzchni adsorbentu oraz wielowarstwowy charak- 
ter adsorpcji rozpuszczalników, Wyprowadzone równania zastosowano do 
opisu doświadczalnych współczynników pojemnosciowych, W oparciu o 
uzyskane wyniki przedyskutowano wpływ wymienionych wyżej czynników 


na mechanizm retencji w adsorpcyjnej chromatografii cieczowej. 


РЕЗЮМЕ 

Зелено уравнение, описывающее зависимость коэффициента 
деления в алсорбционной жидкостной хроматографии от состава би- 
нарной подвижной фазы. Учтено неидеальность подвижной фазы, раз- 
ницы в размерах частиц хроматографированного вещества и раство- 
рителей, энергетическую неоднородность поверхности адсорбента, 
а также многослойный характер адсорбции растворителей. Введен- 
ные ‘уравнения. использовано для описания экспериментальных емко- 
CTHHX коэфрициентов. Опираясь на полученные результаты, продис- 
кутировали влияние вышбупомянутых факторов на механизм ретенции 
в адсорбционной зидкостной хроматографии» 
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